Tunable exciton interactions in optical lattices with polar molecules 
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Rotational excitation of polar molecules trapped in an optical lattice gives rise to rotational exci- 
tons. Here we show that non-linear interactions of such excitons can be controlled by an electric field. 
The exciton-exciton interactions can be tuned to induce exciton pairing, leading to the formation 
of biexcitons. Tunable non-linear interactions between excitons can be used for many applications 
ranging from the controlled preparation of entangled quasiparticles to the study of polaron interac- 
tions and the effects of non-linear interactions on quantum energy transport in molecular aggregates. 



Introduction. The absorption of photons by a solid- 
state crystal gives rise to quasiparticles called excitons. 
There are two limiting models of excitons: Wannier-Mott 
excitons and Frenkel excitons. Wannier-Mott excitons 
occur in crystals with band structure leading to collective 
excitations with an effective radius much greater than 
the lattice constant, while Frenkel excitons are typical for 
molecular crystals, where collective excitations are super- 
positions of elementary excitations localized on different 
lattice sites. These properties lead to important differ- 
ences in non-linear exciton interactions for the two mod- 
els. The interactions between the Wannier-Mott excitons 
are determined by the Coulomb interaction and the phase 
space filling p] , while the interactions of Frenkel excitons 
are determined by shorter range dynamical couplings [2] . 
Multiple experiments demonstrated that Wannier-Mott 
excitons can form two-exciton bound states called biex- 
citons [3HS]- By contrast, despite many theoretical stud- 
ies [MS]) Frenkel biexcitons have eluded the experimen- 
tal observation, with one notable exception |10) . In the 
present work, we show that rotational excitation of ul- 
tracold molecules trapped on an optical lattice gives rise 
to Frenkel excitons with controllable non-linear interac- 
tions. We demonstrate that the exciton-exciton inter- 
actions can be tuned to induce the formation of Frenkel 
biexciton and that the biexciton binding energy can be 
controlled by an external electric field. 

Several experiments have recently demonstrated that 
ultracold molecules can be trapped in a periodic potential 
of an optical lattice [UJ. Such systems can be used for 
the study of quantum energy transfer [TSJ, [T3] , non- linear 
photon-photon interactions |14j , novel quantum memory 
devices |15j and, most notably, for quantum simulation of 
lattice models [TrjH2T] . Although describing different phe- 
nomena, the Hamiltonians presented in these references, 
and in the present work, can be cast in the same form. 
For example, the exciton Hamiltonian discussed here can 
be mapped onto the t-V model, which is the special case 
of the Heisenberg-like models studied in the context of 
ultracold molecules in Refs. [T71[2T]. The key difference 
of the present work from those in Refs. [JTIET] is that we 
explore phenomena associated with the excitation spec- 
trum of the many-body system in the limit of a small 
number of excitations. Because we consider a simpler 
Hamiltonian, our scheme is conceptually simpler, requir- 
ing fewer molecular states and external field parameters. 



We use the rotational states of molecules as a probe of 
the collective interactions, i.e. the experiments proposed 
here can be carried out by measuring site-selective popu- 
lations of the rotational states. This can be achieved by 
applying a gradient of an electric field and detecting res- 
onant transitions from Stark-shifted levels, as described 
inRef. [22]. 

Exciton-exciton interactions in an optical lattice. We 
consider an ensemble of polar diatomic molecules in the 
electronic state trapped on an optical lattice in the 
ro-vibrational ground state. The rotational states of the 
molecules \NMjy) are described by the rotational angular 
momentum N and its projection on the quantization axis 
Mjv- We assume that the molecules are in the Mott- 
insulator phase [TT] and that each lattice site contains 
only one molecule. We consider the rotational excitation 
| JV = 0, M N = 0) \N = 1, M N = 0) of molecules in the 
lattice [23] , For simplicity, we denote the ground state 
of the molecule in site n by \g n ) and the excited state by 
| en). Because the molecules are coupled by the dipole- 
dipole interaction, the rotational excitation gives rise to 
a rotational Frenkel exciton |13j . which is an eigenstate 
of the Hamiltonian: 

JVmol JV mo [ 

£exc = E P„P n + J2 J(n-m)PlP m , (1) 

n— 1 n,m^n 

where J(n - m) = (e n ,g m \V dd (n - m)\g n ,e m ) with 
Vdd(n— m) representing dipole-dipole interaction between 
molecules in sites n and m, Eq is the energy difference 
between the states \g) and |e), and the operators P\ and 
P n are defined by the relations P^\g m ) = S nm \e n ) and 
Pn\e m ) = 8 nm \g n ). The P n and P m operators satisfy 
the bosonic commutation, if n =/= m, and the fermionic 
commutation, if n = m [2] . 

Multiple excitations lead to dynamical exciton-exciton 
interactions described by [2J: 

#dy„=2 E D{n-m)PlPlP n P m (2) 

where D(n - m) = (e n ,e m |V^(n - m)|e n ,e m ) + 

\Vdd(n-m)\e n ,g m ). 
The dipole - dipole interaction operator Vdd(n — m) can 
only couple states of different parity 24J. If \g) and \e) 



2 



are states of well-defined parity, such as the rotational 
states |JVMjv), the matrix elements D(n — m) must be 
zero. 

The inversion symmetry (parity) of molecules on an 
optical lattice can be broken by applying an external dc 
electric field. In an electric field, \g n ) and |e„) are eigen- 
states of the Hamiltonian if" 10 ' = B e N 2 d n ■ £f, where 
d„ is the dipole moment of molecule in site n and £f 
is the electric field vector. They can be expressed as 
l<?) = E N a N\NM N = 0) and |e) = J2 N b N \NM N = 0), 
where and 6jv are determined by the electric field 
strength. 

The Hamiltonian Q can be diagonalized by the 
Fourier transforms: P^(k) — j= — Yl e lkn P^ and 



P{k) 



-ikn 



P n , where k is the wave vector 



of the exciton. This transformation leads to H cxc — 
J2E(k)P^[k)P(k), where E(k) = E + J{k) with J(k) = 

k 

J '{n)e~ lkn , and P'(k) and P(k) create and annihilate 

n 

Frenkel excitons with energies E(k). The interaction ([2| 
in the momentum representation is: 



mol 



2 D(q)P\k 1 +q)P^k 2 -q)P(k 1 )P(k 2 ), 



ki,k2,q 



(3) 

where D(q) = J2 n D ( n ) e ~' qn ■ 

Biexcitons. The exciton-exciton interactions generally 
have little effect on the energy spectrum of two-particle 
continuum states E{k\) + E^fe). However, under cer- 
tain conditions discussed below, non-linear interactions 
may result in the formation of a bound two-exciton com- 
plex, biexciton. The biexciton state is split from the 
two-particle continuum. The splitting is the biexciton 
binding energy. 

Ref. [6] shows that biexcitons can generally form in 
ID and 2D crystals if 



\D\>2\J\ 



(4) 



where D and J are the coupling constants D (n — m) 
and J(n — m) for n — m = 1. In 3D crystals, biexci- 
tons can form if |D| > 6|J| [5]. For ID crystals, the 
biexciton energy in the nearest neighbor approximation 
(NNA) is E b (K) = 2E + D + 4j2 cos ^/ 2 ) [ where K is 
the total wave vector of two interacting excitons and a is 
the lattice constant, and the biexciton binding energy is 
A = (D — 2J) 2 /D. The maximum number of exciton- 
exciton bound states is equal to the dimensionality of the 
crystal, i.e. one biexciton state for ID, two for 2D and 
three for 3D 

For molecules in an optical lattice, the magnitudes of 
J and D depend on the strength of the applied electric 
field and the angle between the field and the molecu- 
lar array (8) (see Figure lc). We calculate these pa- 
rameters for a ID array of 1 S polar molecules (such as 
alkali metal dimers produced in ultracold molecule ex- 
periments) trapped on an optical lattice with the lattice 
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FIG. 1. (Color online) (a) The parameters D and J (in units 
of Vdd = d 2 /a 3 ) as functions of the electric field strength, 
(b) The ratio D/J as a function of the electric field strength. 
The field is perpendicular to the intermolecular axis. For LiCs 
molecules possessing the dipole moment d=5.529 Debye, the 
value £fd/B e = 1 corresponds to £f = 2.12 kV/cm. (c) 
Schematic depiction of the angle 9 between the field (repre- 
sented by blue arrows) and the molecular array (represented 
by red dots), (d) D and J for a ID array of LiCs molecules 
separated by 400 nm as functions of 6 for £/ = 6 kV/cm. 



separation a = 400 nm. Figure lb shows that for a fixed 
angle 9 the ratio D/J increases as the electric field mag- 
nitude increases. For LiCs molecules, the condition Q 
is satisfied for electric fields > 3.6 kV/cm. Note that the 
ratio D/J is independent of 9. 

Frenkel excitons are quasiparticles characterized by an 
effective mass (m e g). The sign of J determines the sign 
of the effective mass [IS] : negative J corresponds to pos- 
itive m c ff and vice versa (see Figure Id). Due to the 
linearity of the Schrodinger equation, a positive poten- 
tial is attractive for particles with negative mass, just 
like a negative potential is attractive for particles with 
positive mass. Because the sign of J and D is the same 
(and consequently the signs of D and m e ff are opposite), 
the dynamical interaction |2]) between excitons in this 
system is attractive. 

To demonstrate the formation of the biexciton and cal- 
culate the biexciton energy, we diagonalize the Hamilto- 
nian H exc + -Hdyn for a one-dimensional array of N mo \ = 
501 LiCs molecules. The matrix of the Hamiltonian is 
evaluated by expanding the biexciton states as 

\* b (K))=T, C * \P\K/2 + k)P\K/2-k)), (5) 

where K = k\ + k 2 and k = [k\ — k 2 )/2 1 and k\ and k 2 
denote the wavevectors of the interacting excitons. The 
Hamiltonian matrix is diagonalized numerically for fixed 
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FIG. 2. (Color online) (a) and (b): Two-excitation spectra 
of a ID array of LiCs molecules on an optical lattice: NNA 
(dashed lines) and exact solutions (solid lines). The shaded 
regions encapsulate the bands of the continuum two-exciton 
states, (c) f)-dependence of the biexciton binding energy A. 
The electric field magnitude is 6.88 kV/cm, 8o = arccos \/2/3, 
0* — arccos \/l/3. (d) Biexciton wave function vs the lattice 
site separation \r\ = \n — m\ of the two excitations for K — 0. 
Inset: Mean width of the biexciton wave function (r) calcu- 
lated as the width of V'jf ( r ) a t half maximum. Numbers at 
each curve indicate the value of D/2.J. 



values of K, which is conserved. Figure 2 shows that 
for 9 = 90° > 9* = arccos(l/v3) the biexciton energy 
is above the two-exciton continuum (binding for particles 
with negative mass), and for 9 = arccos y/2/3 < 9* below 
it (binding for particles with positive mass) . The binding 
energy of the biexciton changes sign at 9 = 9*. It is 
instructive to analyze the biexciton wave function in the 
site representation, which can be obtained analytically 
using NNA: 

mK))= £ e^(" + " l )/ 2 ^(|n-m|)|ptpt i)j 

y/D 2 -4J 2 cos 2 (Ka/2) (2 Jcos(Ka/2) x |rM 



4> K {r) = 



where r = n — m. The biexciton wave function ipxif) is 
plotted in Figure 2d. Figures 1 and 2 thus illustrate that 
the biexciton binding energy and size can be tuned by 
varying the strength and orientation of the electric field. 

Non-optical creation of biexcitons. We now explain 
why it is difficult to observe Frenkel biexcitons in solid- 
state molecular crystals. First, many molecular crys- 
tals, such as anthracene or naphthalene, possess inversion 
symmetry. In these crystals, the constant D as defined in 
the line after Eq.([2]) must vanish and Eq.Q is not satis- 
fied. Second, it is difficult to excite biexciton states opti- 
cally: it was shown in Ref. [7] that the oscillator strength 



for the photon-induced transitions to the biexciton state 
must decrease with the increasing binding energy of the 
biexcitons. Therefore, two-photon excitation can only 
produce unstable weakly bound biexcitons. Third, ex- 
citons in molecular crystals decay via bimolecular anni- 
hilation processes into higher-energy states and subse- 
quent relaxation accompanied by emission of phonons. 
This process is prohibited by conservation of energy in 
an optical lattice with diatomic molecules. Figure 1 
demonstrates that the first obstacle can be removed by 
tuning the electric field. To overcome the second ob- 
stacle, we propose a non-optical method of populating 
deeply bound biexciton states based on the unique struc- 
ture of X S polar molecules. As zero electric field the 
rotational states \g) and |e) are separated by the energy 
Ae e _ ff = 2£? e , while the energy separation between state 
|e) and the next rotationally excited state |/) = \N = 
2,Mn = 0) is equal to Ae e _/ = 4£> e . As the electric 



field increases, Ae e _ s increases faster than Ae e _/. When 
£fd/B e ~ 3.24 (corresponding to Sf ~ 6.88 kV/cm for 
LiCs), Aef-g = 2Ae e - g . At electric fields near this mag- 
nitude, two \g) — ¥ |e) excitons can undergo the transi- 
tion to the |/) state, and, inversely, the \g) —> \f) ex- 
citation can produce a pair of \g) —> |e) excitons or a 
biexciton state depicted in Figure 2. The coupling be- 
tween states \e) and |/) is i? 12 = ^H n ~ m )RnPnPL> 

where M(n - m) = (e„, e m \V d d{n - m)|/„, g m ), and 
the operator R n annihilates the |/) excitation in lat- 
tice site n. The total Hamiltonian describing this three- 
level system is # 9 - e -/ = -ff cxc + H dyn + H 2 + H 12 , 
where H 2 = E f J2 B) n Rn + J2 Jg-f{n — rn)w n Rm and 

Jg-fi n - m ) = (9n,fm\Vdd(n-m)\f n ,g m ). 

In order to calculate the probability of the population 
transfer from state / to the biexciton state, we solve the 
time-dependent Schrodinger equation with the Hamilto- 
nian H g - e -f evaluated in the basis of products of the 
eigenstates of H exc +Ha y n and the eigenstates of H 2 ■ This 
leads to coupled differential equations, which we solve nu- 
merically. Note that the magnitude of J g -f is about ten 
times smaller than J. In the absence of decoherence, the 
\g) — > |/) excitation gives rise to the Frenkel exciton and 
the transition from the / states to the biexciton state 
is a coherent exciton-exciton transition. In the presence 
of decoherence, the exciton states become localized. If 
the decoherence rate is larger than J g ^f/h, but smaller 
than J/h, the \g) — > |/) excitation is localized, while 
the biexciton states remain coherent. Figure 3 presents 
the calculations of the population transfer probabilities 
for both scenarios. The results show that the biexciton 
states can be populated with high efficiency. The equi- 
librium populations (in the limit of large t) depend on 
the relative energies of the / state, the biexciton bound 
state and exciton-exciton continuum states, which can 
be tuned by varying the electric field magnitude. The 
efficiency of the population transfer can be maximized 
if the electric field is detuned far away from resonance 
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FIG. 3. (Color online) Population dynamics for the transition 
from \g) — > \f) exciton (middle panel) and from an / state 
localized on a single molecule (lower panel) to coherent \g) — > 
\e) excitons and biexcitons. The green dashed curves denote 
the population accumulated in the pairs of non-bound \g) — > 
\e) exciton states, the red solid curves the biexciton state and 
the blue dot-dashed curves the / state. The shaded region 
in the upper panel encapsulates the band of the continuum 
two-exciton states. The calculation is for a ID ensemble of 
N mo i = 501 LiCs molecules on an optical lattice with lattice 
separation a = 400 nm. The electric field of magnitude 6.88 
kV/cm is perpendicular to the molecular array. 



when the biexciton population oscillations reach the first 
maximum. Detuning the electric field to low magnitudes 
effectively decouples the / state from the states in the 
{g, e} subspace and interrupts the population dynamics. 
This corresponds to switching off the channel for bimolec- 
ular annihilation of excitons, which is one of the reasons 
of the biexciton population depletion in solids. We have 
confirmed that the calculations with electric fields < 5.0 
kV/cm yield no noticeable population transfer. 

Discussion. We have shown that rotational excitation 
of molecules trapped on an optical lattice gives rise to 
rotational excitons whose interactions can be controlled 



by an external electric field. The exciton-exciton interac- 
tions can be tuned to produce two-exciton bound states. 
A biexciton is an entangled state of two Frenkel excitons. 
The creation of biexcitons as described in the previous 
section and tuning the electric field to the regime of zero 
binding energy can thus be used for the controlled prepa- 
ration of entangled pairs of non-interacting excitons. fn 
order to observe the biexcitons, one could measure cor- 
relations between the populations of the rotationally ex- 
cited states of molecules in different lattice sites using 
the method proposed in Ref. [22] . 

The present work suggests several interesting ques- 
tions. For example, it was recently shown that Frenkel 
excitons in shallow optical lattices can be coupled to lat- 
tice phonons, leading to polarons [55]. Coupling a Frenkel 
biexciton to phonons would produce strongly interacting 
polarons. It would be interesting to explore if these in- 
teractions lead to the formation of bipolarons. 

We have repeated the calculations presented here for a 
system of three excitons and similarly observed the for- 
mation of three-exciton bound states. It would be in- 
teresting to explore the effect of tunable exciton-exciton 
interactions on excitation correlations, both as a function 
of D I J and the density of excitations, to understand fun- 
damental limits of exciton clustering [27 . 

The creation of biexcitons with tunable binding energy 
and measuring quantum energy transport for different 
ratios D/J can be used to study the effects of exciton- 
exciton entanglement on energy transfer in molecular ag- 
gregates j28-3I)|. The ability to tune exciton-exciton in- 
teractions can be used to explore the role of multiple 
excitation correlations on energy transfer in disordered 
systems (the confining lattice potential can be tilted or 
the molecules can be perturbed by a disorder potential 
produced by an inhomogeneous electric field). 
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